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DIFFERENTIAL EQUATIONS OF MOTION FOR COMBINED FIAPWISE 
BENDING, CHORDWISE BENDING, AND TORSION OF 
TWISTED NONUNIFORM ROTOR BLADES 
By John C. Houbolt and George W. Brooks 

SUMMARY 


The differential equations of motion for the lateral and torsional 
deformations of twisted rotating beams are developed for application to 
helicopter rotor and propeller blades. No assumption is made regarding 
the coincidence of the neutral, elastic, and mass axes, and the generality 
is such that previous theories involving various simplifications are con- 
tained as subcases to the theory presented in this paper. 

Special attention is given the terms which are not included in pre- 
vious theories. These terms are largely coupling-type terms associated 
with the centrifugal forces. Methods of solution of the equations of 
motion are indicated, by selected examples. 


INTRODUCTION 


This paper deals with the deformation theory of rotating blades. 

The structural problems of these blades have become more acute in almost 
every phase of aeronautical-engineering application; For example, pro- 
pellers have become larger and thinner, particularly in connection with 
aircraft designed for vertical take-off and landing and short take-off 
and landing, and as a consequence are more susceptible to vibration and 
flutter troubles; helicopter blades are subject to numerous vibration, 
divergence, and flutter problems; and turbine and compressor blades fail 
frequently because of some vibration phenomena. There is therefore much 
interest in the development of a more general deformation theory which 
is fundamental in the structural and dynamic analysis of these problems. 

Although many theories on blade deformation exist, these theories 
either neglect some of the factors of concern or treat them only approxi- 
mately. The factors are many and include variable stiffness and mass 
distributions, noncoincidence of the elastic and mass axes, built-in 
twist, coupling brought about by inertia and centrifugal forces, and so 
forth. In order to give a rough perspective of the scope of available 
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theory, figure 1 has been prepared. The nonrotating -beam cases have also 
been included in figure l(a) . The first two cases refer to uncoupled . 
bending and uncoupled- torsion of- beams without twist, and much work on 
these cases has been done; their treatment is in fact classical. The 
third case is that of coupled bending in two directions and torsion of a 
twisted beam, wherein the elastic and mass axes are noncoincident (see- 
ref. l); the subcase where no twist is present - has also been given exten- 
sive treatment. Additional treatments on general theory of pretwisted 
beams are also given in references 2 to 4. 

The rotating -beam cases which have been treated with the inclusion 
of centrifugal forces are shown in figure 1(b). The first-ncase, that of 
pure flapwise bending of-an untwisted beam, has been considered by several 
investigators; reference 5 is a notable example in which beams of variable 
cross section and different root-end suspensions are analyzed from a 
design-application, standpoint in considerable detail. References 6 to 8 
also give substantial treatments to this case. Pure torsion of untwisted 
rotating beams is treated in references 9 . and 10, and combined bending in 
two directions is treated by various means which include both vector and 
matrix methods in references 11 to l6. 

The case treated in the present paper is shown. in figure 2; this case 
represents the coupled bending in two directions and torsion of a twisted 
rotating beam where the elastic axis, mass axis, and tension axis are not 
necessarily coincident. The tension axis is defined as the spanwise locus 
of the centroids of the cross -sectional area effective in carrying tension. 
The specific purpose of the paper is to develop the differential equation 
of deformation of the blade under the action of various applied loads. 

The development is made along the principles of "engineering” beam theory, 
and secondary effects, such as deformation due to shear, are not - included. 
The theory is therefore intended primarily for blades of moderate to high 
aspect ratios, for which plate bending effects are .probably not significant. 
Special subcases which involve simplifications that appear justified in the 
treatment of certain problems (for example, helicopter blade deformation) 
are then deduced. Finally, methods of solving the equations are discussed, 
and examples are given . 


SYMBOLS 


A 


Bi,B 2 


cross-sectional area of blade effective in carrying axial 
tension 

section constants, see equations ( 9 ) 


ds 


incremental distance measured along beam fiber in undeformed 
state 
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dsi 


incremental distance measured' along team fiber in deformed 
state 


E 

e 


e A 


e 


o 


G 

I l> 1 2 


J 


Young's modulus of elasticity- 

distance between mass and elastic axis, positive when mass 
axis lies ahead of elastic center 

distance between area centroid of tensile member and elastic 
axis, positive for centroid forward 

distance at root between elastic axis and axis about which 
blade is rotating, positive when elastic axis lies ahead 

shear modulus of elasticity- 

bending moments of inertia about major and minor neutral 
axes, respectively- (both pass through centroid of 
cross-sectional area effective in carrying tensions) 

torsional stiffness constant 


k A 




k ml ,k m2 


% 

L z 

M 


polar radius of gyration of cross-sectional area effective 
in carrying tensile stresses about elastic axis 

polar radius of gyration of cross-sectional mass about 
elastic axis = ^nil^ + ^ m 2 ^) 

mass radii of gyration about major neutral axis and about 
an axis perpendicular to chord through the elastic axis, 
respectively 

intensity per unit length of laterally applied aerodynamic 
load in plane of rotation 

intensity per unit length of applied aerodynamic loading 
perpendicular to plane of rotation 

intensity per unit length of applied aerodynamic torque 
loading 


M-i 


resultant cross-sectional moment about major principal axis 


Mg resultant cross-sectional moment about axis perpendicular to 

major principal axis and passing through elastic axis 

M x ,My,M z resultant moments in x-, y-, and z-directions, respectively, 
of M-j_, Mg, and Q moments 
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mass of beam per unit length 

resultant or total loadings per unit length in x-, y-, 
and z-directions 

resultant cross-sectional torque about-elastic axis 
torque loading per unit length 

re suit ant“fc or que loadings per unit length in x-, y-, 
and z-directions 

blade radius 


tension in beam 


, 1 a [* 

J x 


n^mx 


thickness of cross section at any chordwise position 

cross-sectional shears in y- and z-directions 

lateral displacements of beam, in plane of rotation and 
normal to plane, respectively 

coordinate _system which rotates with blade such that 
x-axis falls along initial or undeformed position of 
elastic axis (see figs. 2 and 3) 

blade angle of station x prior to any deformation, positive 
when leading edge is upward 

strain 

m 

strain due to tension, =~- 

ShA 

cross-sectional coordinates; T)-axis lies along major axis, 

| -axis .is perpendicular to major axis and passes 
through elastic axis (see fig. 3) 

values of t] for trailing edge and leading edge of cross 
section 

mass density 

stress 


angle of- twisting deformation, positive when leading edge 
is upward 
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£2 angular velocity of rotation 

o> frequency of vibration, radians /sec 

Primes denote derivatives with respect to x; dots denote derivatives 
with respect to time. 


ANALYSIS 


General Procedure 

As noted from figure 2, the beam is considered as though it were a 
rotating helicopter blade with rotation counterclockwise when viewed 
from above. The x-axis of the x,y,z coordinate axes system shown lies 
outward along the blade and is coincident with the undeformed position 
of the elastic axis. This set of axes moves with the blade around the 
axis of rotation at the given rotational velocity £2, and all deforma- 
tions of the blade are referred to this coordinate system. The blade Is 
considered under the action of distributed aerodynamic loadings in the 
y- and z-directions and under a distributed torque loading about the 
elastic axis, where the Intensities per unit length are denoted by Ly, 

L z , and M, respectively. The tension in the beam Is denoted by T. 

Figures 3(a) and (b) show the coordinate system used for the blade 
cross section and the chosen displacements of this analysis. The q- 
and £ -axes , with the origin at the elastic axis and the q-axis along 
the major axis of the cross section, move with the cross section. The 
blade deformations are denoted by a displacement v of the elastic axis 
In the plane of rotation, positive when in the direction of rotation, a 
displacement w out of the plane ■ of rotation, positive upward, and a 
rotation ( about the elastic axis. The built-in twist 0 and also 0 
are positive when the blade leading edge is up. 

The aim of this analysis Is to derive the differential equation of 
motion In terms of v, w, and 0. The derivation proceeds along the 
following steps: 

(1) The equation for longitudinal strain at any point on the cross 
section is derived In terms of the displacements. 

(2) With the aid of this strain equation the Internal elastic moments 
are derived; these are the resultant moments taken about the q- and £-axes 
and are shown in figure 3(c). 

( 3 ) The transformation is then made of these elastic moments to the 
more easily handled moments which have vectors parallel to the x,y, z axes 
system. (See fig. 3(c).) 
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(V) The equilibrium expressions for these latter moments are- derived, 
and this consideration involves the introduction of the total loading on 
the beam. ■ ■ 

(5) The total loadings, composed of the body forces and applied 
loadings, are derived. 

(6) Steps 3* 4, and 5 are combined to give the final differential 

equations. ■ .... 


Derivation of Strain Equation 

Both longitudinal and shearing strains are, of course, produced 
during bending and torsion of the beam. In general, both_types of strains 
must be considered in determining the resultant forces and moments that 
act on a given cross section. In the treatment— given in the present paper, 
however, an approach is used which requires that— only the longitudinal 
strains have to be considered explicitly. In appendix A a derivation 
akin to an engineering beam theory is given for the strain that "develops 
in any longitudinal fiber of- a twisted beam which undergoes translational 
displacements v and w and a torsional displacement }. The deriva- 
tion is mainly formal, but a physical insight of how strain is developed 
may be gained by considering the possible motions of the imaginary cutting 
planes shown in figure 4 and which are assumed to remain plane during 
deformation. In general, strain may arise from four types of motion: 
pure displacement of the planes toward or away from each other, rotation 
of the planes associated with chordwlse bending, rotation of the planes 
associated with flapwise bending, and rotation of the planes relative to 
one another about the elastic axis to cause beam twisting. With the 
assumption that the cross section is symmetrical about the major princi- 
pal axis, the resulting longitudinal strain is found to be 

e =-£ij, + ( e A " T]) (v"cos (3 + w' f sin [3)+£(v''sin (3 - w"cos (3) + (£2 +T j2_ 

( 1 ) 


The longitudinal stresses follow directly from this equation and are 
a = + ( 'e A - T])(v"cos (3 + w"sin J3) + £(v"sin {3 - w”cos 3) + 

(? 2 + t ) 2 - k A 2 ) p ' 0 ^| 


(2) 
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No further development of .cross-sectional stresses or strains is 
necessary. The consideration of shearing stresses which are associated 
with longitudinal stresses , that is, the shearing stresses which are 
necessary to satisfy equilibrium of an elemental tube, is avoided by 
choosing the elastic axis (defined here as the locus of shear centers) 
as a reference axis; the consideration of shearing stresses due to 
twisting deformation is obviated through use of the assumption that the 
resisting torque includes a St. Venant type torsional term which is the 
same as would develop if the beam were initially untwisted. Both of 
these items will be introduced in the next section. 


Derivation of Internal Elastic Moments 

The stress distributions over the cross section may now be resolved 
into effective internal resisting moments at the elastic-axis position, 
as shown in figures 3 and In order to determine these moments, the 
inclination relative to the elastic axis of the general beam fiber due 
to initial twist and twisting deformation must be considered. The stress 
along this fiber is resolved into two components, one parallel to the 
elastic axis and one in a plane perpendicular to the elastic axis. (See 
fig. 5») From the longitudinal component, the flapwise bending moment Ml 
and the chordwise bending moment Mg are given as follows : 

r\e r t / 2 

M. = - / / <r£ d£ a n (3) 

J \a J - t / 2 


n^le nt/2 

M 2 = - / / crri d£ dr} (k) 

J %e J - t / 2 

where the minus signs have been introduced to make the moments positive 
when they produce compression in the upper and nose fibers, respectively. 

The component in the plane normal to the elastic axis leads to an 
effective torsional resisting moment. The consideration of this compo- 
nent and the addition of the torque associated with St. Venant twisting 
mentioned previously leads to the following equation for total resisting 
torque: 


Q = GJ0* 


f ZS [ cr(|3 + 0)'(n 2 + £ 2 ) d C d n 

J n te J - t / 2 


( 5 ) 
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where a positive torque is associated with a positive 0' . In equa- 
tion (5); the choice of the elastic axis as a reference axis is signifi- 
cant. The elastic axis does not necessarily coincide with the centroidal 
axis, and if the centroidal axis, or for that matter any other axis, were 
used as a reference, then the shearing stresses associated with long!-, 
tudinal' stresses would contribute to the total resisting torque. Such a 
term would have to be included in equation ( 5 ) and would thus lead to 
considerable complication in the analysis. With the choice of the elas- 
tic axis, however, no such term appears, because the elastic axis is 
defined herein as the axis about which the resultant torque of the 
shearing stresses due to longitudinal stresses is zero. 

Since the elastic axis is important, the establishment of this axis 
position will be discussed now. Consider differential longitudinal 
stresses to arise which have a cross-sectional distribution the same as 
that given by equation (2) . The differential longitudinal stresses 

associated with the expression E Jep + e A (V" cos (3 + w"sin p^-k^p'^'J 

are uniform across the cross section and hence do not produce any 
shearing stresses. For the expression 


E jji)(v"c°s p + w"sin p) + + T^p’^’J 

the differential stresses would be symmetrical about the major axis 
(rj-axis) and, since the cross section is assumed symmetrical, would lead 
to a resultant shear directed along the major axis. Only the remaining 
term E£(v"sin p - w”cos p) can lead to shearing stresses which produce 
torque. Thus in order to locate the position of the elastic axis, it Is 
sufficient to consider the beam in bending about the major axis only, 
with a linear stress distribution in the £ -direct Ion, and then to deter- 
mine the shearing stresses over the cross section (in the ip- and ^-directions) 
and the position along the major axis for which these .shearing stresses pro- 
duce- no resultant torque. (See treatments on elastic axis and shear center 
in ref. 17 •) 

The substitution of equation (2) into equations ( 3 ), (4), and ( 5 ) 
gives the following equations for total elastic moments in terms of the 
displacements: 


Mi = EIq(-v"sin p + V’cos p) (6) 

M 2 = El 2 (v M cos p + w"sin p) - Te A - EBgP’^’ (7) 

Q = [gJ + Tk A 2 + EB 1 (p')^' + Tk A 2 p’ - EB 2 p'(v"cos p +-w"sln p) (8) 
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where Iq and I 2 are the principal moments of inertia of the tensile 

carrying area (ig is the moment of inertia about the mean principal axis, 
not about the £-axis) and where Bq and B2 are section constants defined 
by 

B 1 = J le tT - k A 2 ) dT l 

B 2 = / ^ txi^2 + |§ - k/^dT] 

J \e V ' 

In equations ( 9 ) it is interesting to note that if displacements parallel 
and perpendicular to the blade chord had been used Instead of v and w 
the results obtained would differ from those used by other investigations. 
(See appendix A for the results obtained.) Also of interest is the form 
that equation (8) would take if large deflections in 0 were considered, 
but with v = w = 0. This point, which is connected with the possibility 
of a torsional instability, is discussed also in appendix A. 



Moment Transformation 


In the consideration of the equilibrium between moments, shears, 
and tension, it is more convenient to deal with moments that are orien- 
tated parallel to the x-, y-, and z-axes, that is, the moments M x , My, 
and M z shown in figure 3* A simple transformation of the moments Mq, 
Mp, and Q to these newer moments Is thus desired. When the Mq, M 2 ., 

and Q moments are resolved into components in the x-, y-, and 
z-directions, respectively, and use is made of the relations that when 
0 is small 


sin(j3 + 0) = sin 3 + 0 cos (3 

► 

cos(f3 + 0) = cos 3-0 sin 3 
the following relations are obtained: 

M x = Q, + Mq(cos 3-0 sin 3)v’ + Mq(sin 3 + 0 cos 3) v * + 
M 2 (sin 3 + 0 ms 3 )v' - M 2 (cos 3-0 sin 3)w* 

My = Mq( cos 3“0 sin 3 ) + M 2 (sin 3+0 GOS 3) - Qv’ 

M z = -Mq(sin 3+0 cos 3 ) + M 2 (cos 3-0 sin 3 ) + Qw* 


( 10 ) 


(H) 
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Nov viien equations (6), (7)^ and (8) are substituted into equa- 
tions (ll), and all the second-order terms are dropped, the following 
desired equations' for M x , My, and M z in terms of the displacements 
are found: 

M x = [gJ + Tk A 2 + EB 1 (p , )^0' + Tk A 2 f3* - EB 2 p , (v"cos p + w”sin p) - 

Te A v'sin (3 + Te^'cos p (12) 


Eltsin p cos p y" - 

( 13 ) 

M z = ^EI 2 - Elqjsin p_cos p w" + ^EI^ sin 2 p + Elg cos 2 pjv" - 

^Te A + EB 2 p' 0 , ^cos p + Te A 0 sin p (l4) 

Equilibrium Conditions 

The equilibrium of the forces and moments that act on a differential 
beam element is now considered. In this consideration the element is 
formed by slices parallel to the yz-plane, because this choice leads to 
rather simple results. The forces that act on such an element are shown 
in figure (6a); the moments, in figure 6(b). The quantities p x , Jiy, 

p z , q x , qy, and are resultant force and moment loadings, which 

involve both the acceleration body forces and the applied aerodynamic 
loading. The acceleration body forces, due to both centrifugal and 
transverse accelerations, are derived in appendix B. 

Summation of the forces in the x-, y-, and z-directions and summa- 
tion of- the moments about the x-, y-, and z-axes lead to the following 
equilibrium conditions for shear and moment: 


My = ^EI-l cos 2 p +- EI 2 sin 2 p) w M + (eI 2 - 
^Te A + EB 2 P , 0^ sin p - Te A 0 cos p 


T* + p x = 0 
Vy* + Py = 0 ► 



( 15 ) 
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M X T “ VyW* + V z v« + qx = 0 


My 1 - Tw 1 + V z + Qy = 0 


M z 1 - Tv 1 + Vy + q^, = 0 


(l6a) 

( 16 b) 

( 16 c) 


By substitution the shears in these equations may he eliminated to give 
the following basic equilibrium conditions: 




( 17 ) 


M x ' - My'v 1 + M z , w* - q v' + c^v 1 + q x = 0 
M y " - (Tw 1 ) 1 + q.y* - p z = 0 
M z " - (Tv*) 1 + q z I - P y = 0 
Substitution of equations (12), (13)> and (14) into equations (17) gives 

- 1 GJ + Tk A 2 + EB 1 (p , ) 2 ]0' + Tk A 2 {3’ - EB 2 P , (v"cos p + w"sin p)l + 


Te A v"sin p - Te A w"cos p + q y v* - q z w ! - q x = 0 


( 18 ) 


jj^EI A cos 2 p + EI S sin^p)^* + (EI 2 - Eltsin p cos p v ,T - 

^Te A + EB 20 ’ 0 ’) sin p - Tej^ cos pj - (Tw 1 ) 1 + q ^. 1 - p z = 0 (19) 

|^El 2 - EIi^ sin p cos p / + (eIj sin^p + EI 2 cos^p) v" - 

^Te A + EB 2 P , 0 , )cos p + Te A 0 sin - (Tv 1 )' + 5b' - Py = 0 (20) 

where second-order terms have been dropped. Only the loadings now remain 
uo be considered. 
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Resultant .Loadings 

As mentioned-previously, the resultant or total- loadings are com- 
posed of hath, the inertia loading due tor-centrifugal and vibratory accel- 
erations and the applied aerodynamic loading. The. Inertia loadings, 
which consist of the forces and moments which oppose accelerations of the 
blade element's, are derived in appendix B whereas the aerodynamic 
loadings Ly“, L z , and M are left in this symbolic form, since they are 

in the nature of externally applied loadings. 

The procedure that-follows is employed in the derivation of the 
total loadings which is given in appendix B. The acceleration of any 
mass particle on the vibrating, rotating blade is derived and the com- 
ponents in the x-* y-, and z-directions are obtained. These component 
accelerations include terms for the Coriolis force and transverse and 
centrifugal accelerations along with higher order terms. The inertia 
force and moment loadings are then obtained by integrating over the cross 
section; these are simplified by dropping second-order terms and are then 
added with the aerodynamic or applied loadings to yield the following 
desired resultant loadings p x , p^, p z , § x , t3y>_ and q z : 


Px 

= -T' - fi^mx 





% 

- h - m [ v - 

ft 2 (v + e 0 ] 

l] + me|j 

i sin B + fi 2 (cos B - ( sin B)J 


Pz 

=* L z - m(w + 

e^ cos B) 




5x 

= M - + eo)sin 

B + e Q 0 

cos bJ + me(v sin B - ^ cos B) “ 

' (21) 


a2m ( ;k m2 2 " 

Vl 2 ) < sin 

B cos B 

+ 0 cos 2B) - 



= -fi^mex(sin 

B + 0 cos 

P) 



q.z 

= -fi^mex(cos 

B - 0 sin 

P) 




Final Differential Equations 

The substitution of equations (21) into equations (18), ( 19 ) , and (20) 
gives now the desired final differential equations of equilibrium 
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f + 


- |[GJ + Tk A 2 +- EBiCp') 2 ]^ - EB 2 p*(v"cos p + Vein p) 

T e A(v"sin p - w”cos p) + ft 2 mxe(-v , sin p + w’cos p) + sin Pv + 

n2m - Kl) cos 2p + ee Q cos 3 J 0 + n^^'0 - me(v sin jB - w cos (3) = 

M + (lk A 2 p')' - S&l (k^ 2 - ]%,-[ 2 jsin p cos p + ee Q sin p^j (22 


j^EIp cos 2 p + Elg sin 2 p)w" + (eI 2 - EI^J sin p cos p V* - Te A $ cos P - 
EBgP'^'sin pj" - (Tw 1 ) * - (ft 2 mxe0 cos p)’ + m(w + e$ cos p) = 

L z + (Te A sin P) " + (fAnxe sin p) 1 (23) 


|~(eI 2 - EIp) sin p cos p w" + (eIj_ sin 2 P + EI 2 cos 2 p^v" + Te A 0 sin p - 
EB2P t 0 t cos pj " - (TV* ) 1 + (fl 2 mxe(/f sin p) 1 + sin |3 + m(v - e0 sin p) - Q?i 

L y + ( Te A cos p)" + (ft 2 mxe cos p J 1 + A(e 0 +e cos p) ' (24) 


In these equations, the terms that are not included In previous 
theories have been underlined, (it should be noted that ref. 18 suggests 

the existence of the centrifugal coupling terms ft 2 mxew’ and ( ft 2 mxe0 ) 
for the case of a blade with p = 0.) Many of these new terms are 
coupling-type terms which are associated with the centrifugal forces. A 
number of the new terms take the forms of lateral and twisting loadings 
which tend to deflect the blade even in the absence of externally applied 
loads; because of this fact they have been written on the right-hand 
sides along with the applied loads Ly, L z , and M. 

Importance of coupling forms .- It is beyond the scope of this paper 
to make a complete evaluation of the magnitude of each of the coupling 
terms and their influence on such phenomena as, say, blade vibration. A 
rough idea of their importance can be gathered by analyzing each term in 
detail and comparing them to other terms of known importance. As an 
example, consider (Te A 0 cos p)” in equation (23). If, for simplicity 
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in the consideration, e^ and m are considered independent of x, and 
if (3 = 0, then this term would become oh expansion 


-ft 2 me A 0 + 2T'(e A 0)» + T(e A 0 } " 

These terms may now be compared with other terms of known importance. 

Thus, if e A 0 is interpreted as a deflection, then the last two terms 
of this expansion are seen to be exactly comparable to the expansion of 
the well-known and important term (Tw 1 ) 1 = T’w‘ + Tv” . Also the first 

term -fl^m e A 0 is directly comparable with the coupling term of estab- 
lished importance in the consideration of coupled bending-torsion of 

vibrating beams, namely, o>2me0. A similar consideration may be given 
each of the other new terms in the final equations. The examples to be 
presented later- will also give an indication of the importance of some 
of the terms . 

Boundary conditions .- Since the problem is, of course, a boundary- 
value problem, some mention of the associated boundary conditions should 
be made. The most general case involves ten boundary conditions - two 
conditions associated with equation (22), and four each. for equations (23) 
and (24). These conditions all evolve from the type of— constraint^ that - 
is imposed on the ends of the beam on the displacements v, w, and (A . 
the moments My and M z (eqs. (13) and (l4)), the torque Q (eq. (8;), 

and the shears Vy and V z which are defined by equations (l6b) and (l6c) 
For a completely fixed end, the boundary conditions would be 

v = w = 0 = v' = w' =0 
For a free end, the conditions would be 

Q = My = M z = Vy = V z = 0 

In this case, the alternate choice of Mp = M 2 = 0 (see eqs. (6) and (7) 
may be used instead of My = M z = 0. 

For the present - system, an end cannot be-referred to simply as pin- 
ended; care must be taken to specify the number and directions of pins 
present. For example, an end may have a pin which runs parallel to the 
y-axis but may be fixed as regards displacements v and 0. In this 
example, the moment My has been relaxed and the boundary conditions 

become 


v = w = v * = 0 = My = 0 


Other pin-ended conditions follow in a similar fashion. 
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Elastically restrained ends may be handled in a manner similar to 
that used in the usual treatment of beams , except that the direction of 
the restraint must be properly taken into account. 

Subcases of the general theory .- In actual applications, several 
special subcases of the present theory are frequently encountered. Three 
of these subcases are described as follows: 

Case I: The built-in twist is zero, p = 0. 

Case II: Bending in two directions without torsion, 0=0. 

Case III: Torsion only, v = w = 0. 

The equations of motion which result when the theory Is reduced to these 
special cases are as follows : 

For case I, £ = 0, 


-[(GJ + Tk A 2 )0' - Te^w " 4* ft 2 mxew I + ft 2 m ^ k^? 2 - 2 + ee o ^0 + ml^ 2 ^* mew= 


(eIiw m - Te^)" - (Tw 1 ) 1 - (n^nxe0)' + m(w + e$) - L 2 


(e^v") ” - (Tv*) 1 - ft^v + mv = ly - (Te^) 1 * + (si 2 mxe) t + + e) 


For case II, 0=0, 


|j(EIi cos^p + EI 2 sin^p)w" + (EI 2 - Eli) sin p cos p v"J " - (Tw r ) f + mw = 
L z + (^ e A sin P)" + (Axe sin p) f ( 28 ) 

£(eI 2 - Eli) sin 3 cos p w" + ^ET-l sin^p+EIg cos 2 p)v'^J 11 - (Tv 1 ) 1 - ft^mv-h mv = 
Ly+ (Te^ cos p) n + (ft^mxe cos p) 1 + ft^m(e 0 + e cos p) (29) 
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For case III," v = w = 0, 


[GJ + Tk A 2 + EB 1 (j3 , ) 2 ]0 , | + fAnj^k^ 2 - k^^cos 2p + ee Q cos pj0+mk m 2 0 = 

M + (Tk A 2 j5 ' ) 1 - - k^^sln p cos p+ee Q sin pj ( 30 ) 


Associated Energy Equation 


As a check on the derivation presented herein, the differential equa- 
tions of equilibrium were derived by a completely different approach that— 
involves energy principles. In this energy approach the strain energy of 
deformation, in which the- longitudinal strains were arrived at in a dif- 
ferent manner than that— presented herein, and the work performed by all 
the forces present were considered. The equation for total potential 
energy of the system is, of course, closely allied to the differential 
equations and is of intrinsic value in the treatment of the problem of 
the present paper from an energy standpoint. For completeness, the equa- 
tion is presented here without derivation. If U denotes the total strain 
energy and V represents, the work performed by the centrifugal body forces 
and the applied loading, then 


U - V = 


i J |EI 1 (v"sin p - w"cos p) 2 + El 2 (v"cos p + w"sin p) 2 + 
|gJ + EB-lCp 1 ) 2 ] (0’ ) 2 - 2EB 2 (v”cos 3 + V’sin 3)3’ 0‘| dx - 
J 0 + ( v ’) 2 ] + e A (v n cos P + w"sin p) - 

e A 0(v"sin p - w J, cos p) - ) 2 + p^'J ► + 

xejj(v’cos p + w’sin p) + 0(v’sin p - w'cos p) ] + + 


n 2 m I 

(e cos p + eo - e0 sin p)v + 
ee 0 sin p 
Pyv + p z w + q0)dx 


[-( 


k m2 2 - k ml 2 ) sln ? cos P ’ 


*]$ " “ k ml 2 ) cos 2 P + ee o cos J^0 2 - 

\ 


( 31 ) 
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where Py, p z , and q are the applied loadings considered in a static 
sense. Application of the minimum variational principles to this expres- 
sion yields differential equations which are the same as those given by 
equations (22), (23), and (24) when 0, v, and w are not time dependent. 
The variation also gives the boundary conditions on moments, shears, and 
torques that are for stipulated geometrical conditions at the ends of the 
beam and are given in the section entitled "Boundary Conditions" . The 
differential equations applicable when time dependent motion is considered 
may now be obtained by writing the loadings p^., p z , and q as the sum 

of the applied loadings and vibratory inertia loads (as obtained by 
D 1 Alembert ’ s principle ) ; specifically, 


Py = Iy - m(v - e<j) sin (3) 
p z = L z - m(w + e0 cos J3) 


q = M + me(v sin p - w cos p) 



(32) 


The substitution of these loadings into the differential equations obtained 
by the variational process then gives equations identical to equations (22) 
( 23 ), and (24). 


METHOD OF SOLUTION AND EXAMPLES 


The general differential equations of this paper cannot be solved 
exactly, and it is therefore necessary to resort to some approximate means. 
Two means, which differ in manipulations but which yield equivalent results 
will be considered here - one is a modified Galerkin type procedure, the 
other a Rayleigh-Ritz procedure. 

In the modified Galerkin procedure, the displacements are assumed to 
be expressed In finite series of known functions as follows: 


0 = a^^ + a 2 0 2 + 
w = b^w-^ + b 2 w 2 + 

v = c l7l + C 2 V 2 + 



• • ' b q w q 




(33) 


where a n , b n , and c n are unknown coefficients, which may be time 
dependent, and 0^, w n , and v n are chosen modal functions which satisfy 



18 


mCA TN 3905 


the geometrical boundary conditions of the problem. A set of linear 
simultaneous equations which allow the evaluation -.of the unknown coeffi- 
cients is now formed through application of the following equations: 


T 0 n D 1 (^,w,v)dx + Q0 n l = 0 

J 0 -*0 


(n = 1; 2 } . • • p) 


pE R _R 

/ v n D 2 (0,w,v)dx + M y w n ' + V z w n = 0 
Jo 0 0 


nR -| R .. -jR 

J v n D 5 (0,w,v)dx + M z v n 'J ^ + VyvJ^ = 0 


(n = 1, 2, 3, 


(n = 1, 2, 3, 


(34a) 


. . . q) 

(34b) 


. . . r) 

(54c) 


where Dj_(0,w,v), Dg(0,w,v), and D^(0,w,v) denote, respectively, equa 

tions (22), (23 ) y and (24) with all the right-hand" terms transferred to 
the left-hand side; Q, My, and M z are given by equations (8), (13) , 

and (l4) ; and V z and Vy are found from equations (l6b) and (l6c). 

For static or other externally imposed loading conditions, the set of 
equations formed by equations (34) will be nonhamogeneous ; the solution 
for the unknown coefficients and, henceT-the displacements then proceeds 
by ordinary means. For a characteristic -value problem, such as free 
vibrations, the set— of equations will be homogeneous, and the vanishing 
determinant will yield the characteristic frequencies. 


This procedure is called a modified Galerkin procedure because of 
the presence of the boundary terms in equations (34) . In general, it 
would be very difficult ter choose modal functions in equations (33) so 
that both the geometric and "force" (shear, moment, torque) boundary con- 
ditions are satisfied, which is necessary in the usual Galerkin approach. 
Therefore, the force boundary terms are added in equations (34), and, 
thus, the stipulations on the chosen modal functions are relaxed so that 
only the geometric conditions need be satisfied. It should be mentioned 
that in actual applications these force boundary conditions do not- have 
to be evaluated, since they can always be cancelled by integrating certain 
of the terms of- the integrals by parts. A justification of this procedure 
can be obtained from the Rayleigh-Ritz procedure. 

In this Rayleigh-Ritz approach the displacements are- expanded in 
series as_ before (see eqs. (33) )> with the same requirement on the geo- 
metric boundary conditions. These expansions are substituted into the 
energy equation ( 31 ) , and a minimization is then made with respect - to the 
unknown coefficients; that is, 
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d(U - v) _ 
^ a n 


(n = 1, 2, 


p)l 


d(U - v) = 

Sbn 

d(U - v) = 

Sc n 


(n = 1, 2 , 
(n = 1, 2 , 


9.) ^ 


(35) 



The resulting equations will he Identical with those of equations (3*0 . 

As a matter of fact, equations (35) may "be used to prove the validity of 
equations (3*0 (through appropriate integration by parts) and thus to 
show that the two procedures outlined here are really equivalent. Because 
of the confusion that often arises about boundary conditions when the 
Galerkin process is applied, the safest procedure is to use the Rayleigh- 
Ritz approach. 

Two examples are now given to show the application of the two pro- 
cedures discussed and to show further the importance of some of the 
coupling terms. 


Example 1 

In example 1, it is desired to determine the natural frequencies of 
a rotating cantilever blade having zero initial twist or blade angle of 
attack (p = 0) and with v = 0. The Rayleigh-Ritz procedure will be used. 
Assume that the displacements 0 and w are given as follows: 


0 — ap0p + a 202 t- a ^0^ 

w = b^_w^ + b£W2 + b^Wj 

where 0 n and w n are the natural uncoupled modes of the nonrotating 
beam, in torsion and bending, respectively, and are expressed in terms of 
a -unit tip displacement.' Equations ( 36 ) are substituted Into equation (31) 
with v = p = 0 , and a minimization is made with respect to and b n , 

and then use is made of equations ( 32 ), with the applied loadings 
Ly = L z = M = 0. The resulting equations are then reduced as follows. 
First, the relations which apply to uncoupled vibrations are used, namely 


(36) 
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mw m v n ^ = 0 


= M 


■n 


(m ( n) 
(m = n) 


0 


dx = 0 


= / raw ^dx 

^ Jo n 


= “n 2 ^ 


mk A 2 0 m 0n dx = 0 


= N. 


n 


(m i n) 
(m = n) 


> ( 37 ) 


(m / n) 
(m = n) 


^ = 0 


p R 


- ^m^n 2 (m “ n > 


(m ^ n) 


= t% 2 H n 


J 


where and p n are the natural-circular frequencies of the uncoupled 

bending and torsion modes, respectively. Second, the constant terms are 
dropped, since they apply to static deformation of the blade. And third, 
motion representative of simple harmonic motion is considered by setting 


a n = a n s ^- n ttfk 
b n = b n sin tut 

The final equations will be homogeneous, and the- vanishing of the deter- 
minant of the a n and b n coefficients will define the natural coupled 
rotating frequencies. This determinant is 
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(“l 2 “ <n2 ) M l +A li 

A 12 

*13 

-B-) 1 7 

-Bj>7 +Q^Dpi 

-B ^ jj-cc^C 51 

A 12 

(cD 3 2 -gl^) tfe+Ago 

*23 

— B p-<D^C-^ p-f ft^D-] p 

— B pp-fB^C pp+Q^Dpp 

-Bjg-OJ^C Jg+fl^bjg 

a 15 

*23 

( m 3 2 ' tc 2 ) M 3 +A 33 


-B 25 -<n?C 2 yMi 2 D 23 

^ aj j -a^C j 

— B-| i 1 +n 2 P 1 ■) 

-Bt p+Q^St p 


Ei-fEn-t-Q^’ n 


e 1J +q ^13 

-B ppKD*~C 

“B 


Big* Q %to 

(tig 2 '<o 2 )yp+E^j+Q^F g j 

E25+0% 2 j 

-B ^-^d^C^X+Q^) 

-B 

-B33-aft33+Q 2 D J 3 

Ei^+Q 2 ?^ 

^23+^23 

(u 3 a ^ 2 )K 3 +E 33 +n ^ 3 3 


where the coefficients are 



No numerical evaluation of this example Is made, hut the Importance 
of some of the new coupling terms of the present theory Is realized hy 
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comparing the associated coupling terms which appear in the determinant . 
For example , ‘ fl 2 P -| 1 , which is new, may be compared with ai 2 C]_]_, a con- 
ventional term of established importance. The frequencies ft 2 and cd 2 
may be of the same order of magnitude; the factors under the integral 
are the same except for w^_ and xwp 1 . These may be shown to be of the 
same order of magnitude, however, and so both coupling terms are com- 
parable in magnitude. This is another indication, therefore, that the 
new coupling terms are of significance in the treatment of rotating 
blades. - - 


Example 2 


The case treated for example 2 is a uniform blade with an offset 
mass at the tip and with p, v, e A , and e Q ' equal 0. (See top sketch 
in fig. 7-) This system may be regarded as having a behavior similar lo- 
an actual helicopter blade which is heavily loaded along .the leading 
edges, as depicted in the lower sketch of figure J. An approximation to 
the first _ two natural frequencies will be found by the modified Galerkin 
procedure. The differential equations for shear, moment, torque, and ten- 
sion (eqs. (25), (2 6), (l6), (13) , and (8) with p = v = e A = e Q = 0) 

applicable here are 


- |^GJ + Tk A ^0’J + U^mxew* + fi'un - k^ ^0 - o^mew 4- mk^^ = 0 

(EIpw n ) - (Tw 1 ) 1 - (fl^mxe^)' - cn 2 (mw + me0) = 0 


( 39 ) 


My = EX 1 w" 

V z = -(e^w")* 4 — Tw* + G 2 mex0 
Q = (GJ + Tk A 2 )9S* 
p R 

T - fimnx dx + 

J x 


(4o) 


In these equations the mass per unit length m should be interpreted to 
apply to both the distributed mass of the beam and to the concentrated 
offset mass. This concentrated mass can be handled, for example, by 
imagining the mass to be distributed over an infinitesimal length, say 7\, 
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with intensity M/A, then to perform the integrations of the Galkerin 
process as though all. quantities involved in the integrations are uniform 
in this interval, and finally to squeeze the interval down to zero length. 
An equivalent procedure is to take the concentrated mass into account 
through means of the Dirac delta function. 

With this in mind, equations (39) and (40) are substituted into 
equations (34a) and (34b); the terms 0 n (GJ + Tk A 2 ^0 T and 

w n £(Tw‘ )• - (n 2 mxe0) f J are integrated once by parts, the term w n (Elqw") " 

is integrated twice by parts, the boundary terms cancel, and the results 
are the following two equations : 


R 


r r - 

(GJ+Tk A 2 ) 0 , 0 EL I dx+ J ft 2 mxew , +Q 2 m^k^ 2 -k^ 2 ^0-Oi 2 ^mew+mjl% 1 2 0j 


0 n dx = 0 

(41) 


R R ] 

J EIqW n " dx + J ^Tw* + fl 2 mxe0^ w n * dx - a> 2 


(mw ■** me0)w n dx = 0 
(42) 


Approximations to the displacements are now chosen as 


0 = aq0q 

w = b fl 


(43) 


where 0q and wq are the fundamental uncoupled vibration modes of the 
beam without the tip mass, in torsion and bending, respectively, and each 
is given in terms of a unit tip displacement; these functions satisfy the 
equations 





mk m 2 0 1 2 dx = [i-, 2 ^ 2 | 




mwq 2 dx -= 60 q 2 




2 . \ 
2 ml ^ 2 R 2 j 


( 44 ) 



J 



24 


MCA TN 3905 


where m refers only to the mass per unit length of the beam. Equa- 
tions (43) n.ov substituted into, equations . (4l) and (42), use is made 
of equations (44) and the equation for tension (applicable here) 

T = i sAn^R 2 - x 2 ) + and the integrations are performed with due 

cognizance being given the offset tip mass; the values found for the 
integrals are as follows : _ 



GJ^dx = p 1 2 mk m 2 | 
U 2 mxew 1 '0 1 dx = l.^fl^a 
fl 2 m (h^p 2 - 2 j 0i 2 dx = U 2 Ma 2 
mew-^j_ dx = Ma 
ml% 1 2 ^]_ 2 dx = mk^ 2 — + Ma 2 
EIi(wi") 2 dx = oo-j 2 ^ 

T^Wi') 2 dx * a 2 ^(mR + | m) 
mw 1 2 dx = 4 - M 


(45) 


(The contribution of Tk^ 2 is found neglible in the first' integral.) 
The final equations are given now in nondime ns ional form 


M-l 2 mRkm 2 
ojq 2 ' 2Ma 2 


+ A 2 ' cd 2 

cop 2 aq 2 \2Ma 2 


+ 1 


+ 
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-Y 




inR + fl 2 At- + LmR | cu^ A + mR\ ^1 _ 

4m 0^2^ 15M/ “ mi 2 \ W a 


The vanishing of the determinant of these equations will yield the 
natural frequencies. Specific evaluation has been' made for the fol- 
lowing values: 



= 1/0T02 



0.152 


and the results are shown in figure J, Results are shown for three con- 
ditions: (a) for vibration of the beam with offset mass in pure bending 

or pure torsion, (b) for vibration in the conventional coupled sense but 
with no centrifugal-force coupling (that is, the underlined terms of the 
equations of this example are emitted), and (c) for vibration with all 
coupling terms considered (the curves labeled complete in fig. 7 )* It 
may be seen that the inclusion of centrifugal-f orce coupling, which is 
new in this paper, has a very pronounced influence on the vibration 
characteristics. In fact, this centrifugal-force coupling seems to be 
of the same general strength as the conventional offset mass coupling 
normally considered. 


CONCLUDING REMARKS 


The differential equations of deformation of a twisted rotating 
blade under arbitrary loading have been derived. A special feature of 
the analysis is the consideration of the noncoincidence of the tension, 
mass, and elastic axes, whereby coupling is allowed to occur in a more 
general form than heretofore considered. This theory can therefore be 
used to solve problems in which these effects are significant and to 
evaluate the less complete theories which perhaps are more easily applied. 

The aerodynamic loading considered was left in general form. It 
can be replaced, however, by whatever terms are appropriate for the case 
being treated. If these aerodynamic forces are taken as zero, then the 
loadings for free vibration result. 
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Because of its usefulness, the energy equation for blade deformation 
is also presented. Two methods of solution are discussed: one is a modi- 

fied Galerkin process 'which makes use of the differential equations; the 
other is a Rayleigh -Ritz procedure which makes use of the energy expres- 
sion. Two examples are presented which illustrate the application of 
both of these procedures. 

A discussion of the importance of the new coupling terms that arise 
is made and their significance is further brought out by the examples . 

The indications are that the centrifugal-force coupling may be as impor- 
tantras the mass coupling that is normally considered in beam analysis 
and therefore should be included as a regular part in the treatment - of- 
rotating blades. ' r 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., October 5, 1956. 
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APPENDIX A 

DERIVATION OF LONGITUDINAL STRAINS 


In this appendix the equation Is derived for the longitudinal strain 
that is developed "when a twisted beam undergoes both lateral and twisting 
deformations. Consider an imaginary plane to cut through the beam per- 
pendicular to the elastic axis; the location of a fiber f of the cross 
section, both before and after the deformation, may then be given according 
to the following sketches: 



The x-axis is normal to the paper and is made coincident with the unde- 
formed position of the elastic axis. In terms of the distance rj along 
the major axis, the distance £ perpendicular to this axis, and the 
built-in twist p, the initial y- and z-positions of the fiber and the 
rate of change of these positions with respect to the x are 


y = T] cos p - £ sin (3 
z = r] sin p + £ cos P 


y 1 = --qP’sin p - £p’cos p = -P T z 
z ! = T}P 1 cos p - £p’sin p = p*y 

Now consider beam displacements to occur so that the point of Intersec- 
tion of the elastic axis and the cutting plane moves the distances u, 
v, and w in the directions of x, y, and z, respectively, and so that 
the cutting plane remains perpendicular to the elastic axis and rotates 
around it by an angle 0 . Then if the usual small -angle assumptions 
that cos(p + 0 ) = cos P - 0 sin p and sin(p + 0 ) = sin p + 0 cos p 
are made, the new position of the fiber is defined by the following 
equations : 
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Xq = x + u - v'^ 

y 1 = v +- T}(cos p - 

Zj_ = w + T)(sin p + 

and the derivatives 
as follows: 


- vj - w' ^z^ — w^=x+(i — v'y - w'z 
0 sin p) - £(sin p + 0 cos p) = v + y - z0 . 
0 cos p) +• £(cos p - 0 sin p) = w + z +• y0 


(A3) 


J 

of x^, y^, and zq with respect to x are given 


x-^' = 1+ u’ - v"y + v’p’z - w M z - w'p'y = 1+ u* - y(v" + p T w f ) - z(w"- P'v')! 


y^ 1 = v ’ - p , z-p , y0-z0 , = v , -yP’0-z(p f + 0 r ) 


UA4) 


z 1 , = w’ + P , y-P , z0+y0* = w* + y(p* + 0') - zp’0 


The longitudinal strain that is developed in a fiber may he found 
from these equations "by considering the amount an elemental fiber of 
length ds changes in length as a result of deformation. In terms of- 
the differential components of length in the x-, y-, and z-directions, 
the final length dsq of a fiber is given by the following equation: 


dsj_ 2 = dxj_ 2 + dy-j_ 2 + dzj_ 2 


Thus 


(Irf = ( x i ') 2 + ( n ') 2 + ( z i ') : 


(A3) 


(A6) 


which becomes with the use of equations (A4) 
2 


(ix^) = I + 2 - y(v" + P’w') - z(w” - P'v* )j +- z 2 p' 2 - 

2zP , (v* - yp'0 - z0') +-y 2 p ,2 + EyP'Cw' + y0 f - zp*0) (A7) 


or 

dsi 

dx 


1 + (y 2 + z 2 )(p’ ) 2 + 2ju' - yv" - zw 1 ’ + (y 2 + z 2 )p , 0 , j| 


1/2 


(AB) 
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Nov the analogous equation for the original length ds may be determined 
directly from equation (A8) by letting u=v = v = 0= O. Thus 


[l + (y2 + z 2)(p.)2] l/2 


(A9) 


The tensile strain in the fiber can nov be written 


e - 


dsj_ - ds 
ds 


ds- 


ds 


- 1 


= U + 


— ][u* - yv" - zv" + (y 2 + z 2 )p'0*jj- 


1/2 


|_1 + (y 2 + z 2 )(p' 

1 1[ u . - yv” - zv" + (y 2 + z 2 )p , 0’j 


1 + (y 2 + z 2 ) P' 2 


- 1 


(A10) 


Nov with ( y 2 + z 2 )(p 1 ) 2 « 1 (say less than 0.03 to 0.04, vhich is 
generally the case) ; then 


e = u 1 


- yv" - zv" + (y 2 4- z 2 )p'0’ 


(All) 


Use is nov made of equations (Al) in order to express the strain in 
terms of the cross-sectional coordinates rj and £ as follovs: 

e = u* - r] ( v" cos p + v"sin p) - £(-v"sin p + v" cos p) + ^r| 2 + £ 2 j p 1 0 1 

(A12) 

It is convenient nov to eliminate the strain component u'; this 
is done by making use of the equilibrium condition that the integral of 
the longitudinal stress over the cross section must be equal to the total 
tension. Thus, with the assumption that the cross section is symmetrical 
about the major principle axis, the folloving equation applies: 

pt / 2 

T = E / / e d£ dq 


■t/2 


P t / 2 

= E / / lu r - rj(v"cos p + v"sin p) - £(-v"sin p + v"cos p) + 

Ate J -*/2 

(n 2 + £ 2 )p'(zf']a£ ai) (A 13 ) 
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or .. : .. 

T - E Jau* - Ae A (v"cos J3 + w"sin p) + Ak^P'^H (Al4) 

This equat ion yields 

u* = + e A (v"cos p + w 11 sin p) + k^P'^ 1 (A15) 

If T/EA is denoted "by the tensile strain eq, then equation (A15) com- 
bined -with equation (A12) yields 

e = + (e A *- T])(v n cos p + w"sin p) + £(v"sin p - w’^os p) + 

(tj 2 + g 2 -- l^ 2 )p»0' (A16) 

which is the complete expression for the strain of— any fiber in the 
cross section. 


Expression for Strain in Terms of Displacements in 
and Normal to the Blade Chord 

Because of apparent anomalies that have arisen in the pastq it is 
of interest~t'o show the development of the strain in terms of displace- 
ments in and normal to the blade chord. These displacements have been 
used by some investigators and are shown in relation to the displacements 
of this paper in the following sketch: _ 



These displacements are related according to the equations 


v = vq cos P - wq sin pi 


w = Vq sin p + wq cos p 


(A17) 
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If two derivatives of these equations are taken, and due care is given 
the fact that (3 may vary with x, then the following combination of 
these derivatives may he found: 


v" cos p + w" sin p = vp" - 2(3 , w 1 I - fS"wp - (p , ) 2 v 1 

p 

-v"sin p 4- w"cos p = wp" + ^'vp 1 + p n Vj_ - (p*) w-^ 


(ai8) 


j 


If these expressions are now substituted in equation (Al6), then the equa- 
tion for strain becomes 


e = + 





+ 2P'v 1 l + P'V^ - 


(0 r ) 2 ^ ± 


+ (rfi + £ 2 - k A 2 ^p ! 0 l 


(A19) 


The interesting fact to be noted here with respect to vp and wp is 
that the strain is not proportional only to the derivatives V]_" and wp" , 

as might first be supposed, but rather to additional terms which involve 
the rate of change of initial twist . 


This observation has some significance. If the analysis of this 
paper had been carried through in terms of displacements vp and w-|_, 
then the moments and all the results would be in terms of the group of 
terms appearing on the right-hand side of equations (Al8), and the final 
result would be in agreement with the results presented in reference 11. 
This consideration indicates that the treatment in reference 1 may be 
subject to correction, since the assumption is made in this reference 
that the moments are proportional to vp" and wp" . Evidently this^ 

assumption is not valid when the beam has a finite rate of change of 
Initial twist. 


Strain When 0 Is Not Small 


The preceding derivation was made on the assumption that the dis- 
placements u, v, w, and 0 were small, and a linear equation for 
strain results. An interesting observation may be made, however, for 
the case of v = w = 0, but with large displacements in such that 
nonlinear effects are present. Thus, if the derivation is repeated with 
<f considered large, and with v = w = 0, then equation (All) would take 
the form 


e 


= u 1 + ^y 2 + 
= u r + ( q 2 + 


Z 2 ) jp'0’ + |(0') 2 

e 2 )[p , 0' + 1(0') 2 ] 


(A20) 
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The strain component u' is now handled as before j and the -equation for 
strain becomes - 


e 





(A21) 


This equation may be used to demonstrate the fact pointed out - in 
reference 19 that - pretwisted blades may pos~sess a torsional instability. 
Thus if use is made of equation (5); and the stress is based on equa- 
tion ( A21 ) , then the following equation for internal elastic resisting 
torque results : 


0 = GJ0 1 

+ Ik A a (0' + 0’) 

+- EB]_ 


= ‘G J0 ' 

+ Tk A 2 (fs’ + 0') 

+ EBl 
2 

(P* 


+ |(0') 2 J(P' +-0') 

+ 0') 3 - (P') 2 (P' + 0')] 


(A22) 


Observations similar to "those made in reference 19 may be made about 
equation (A22) with respect to the possibility of torsional instability; 
as a matter of fact, if the - net tension over the cross section is zero, 
the equation becomes similar to equation ( 2 ) presented and discussed in 
reference 19 . 
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APPENDIX B 


DERIVATION OP ACCELERATIONS AND INERTIA LOADINGS 


The purpose of this appendix is to derive the acceleration of a 
mass particle of the twisted rotating blade and to derive the inertia 
loadings from this acceleration. In order to do this it is convenient 
to introduce the fixed X,Y,Z coordinate axes system shown in the fol- 
lowing sketch: 



This sketch also shows the rotating x,y,z axes system (used in the body 
of the report and in appendix A) and the deformed positions x^, y^, and 

z-j_ of the mass particle (see eqs. (Aj)). The azimuth position of the 
rotating system relative to the fixed-axes system is denoted by Dt. 

In terms of the unit vectors i, j, and k, the vector r may be 
written 

r = ^xj_ cos ftt - y 2 sin Dtj i + ^x^ sin Sit + y 2 cos Dtjj + z^_k (Bl) 


Differentiating with respect to time gives the velocity vector of the 
mass particle 



cos ftt - ftx^ sin fit - yg sin ftt - fiy2 cos 
sin Dt + fiX]_ cos ftt + yg cos Sit - Dy 2 sin 


ftt^i + 

fttjj + z-jk 


(B2) 
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and differentiating once again gives the acceleration 
r = (x x - ay 2 ) cos at - a^xq - ay 2 ) sin at - (y 2 + ax^ sin at - 

^x’i - ay 2 ^sin at-+ a ^X]_ - ay 2 ) cos at + 
($2 + cos at-= a(y 2 +■ ax L J sin at) j + zqk 


i(y 2 -h fix^cos atj 


i + 


(B3) 


The components a x , a y, and a z of the acceleration vector in the 
X-, y-, and z-directions may be-- found from this equation by letting 
t = 0, thus 


r l = (*1 " ^ 2x i - 2ay 2 )i + ^y 2 - a 2 y 2 +--2ax^ J + ^k 


= a x i + a y j + a z k 


(B4-) 


Use is now made of equations (A3) and the fact that y 2 = y± + e Q . 
Thus, 


X]_ = x + u - v'y - w'z 


y 2 = v + y - z0 + e c 


^1 = w + z + y0 


The first two derivatives with respect to time are given as follows : 


* ■ * • • . ^ 

= u - v'y - w'z 

y 2 - v _ 2 0 

Z]_ = W + y$ 


(B5) 
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x-L = u- v I y-w , z 

y 2 = V - z0 > 

•• • • V 

zq = w + y0 


(B6) 


Substitution of these equations into equation (B 4 ) yields the desired 
acceleration vector , with the following components in the x-, y-, and 
z -directions : 


a x - u - v’y - w’z - ft 2 (x + u - v’y - w’z) - 2fl(v - z$) 


= v - 20 - Q 2 (v + y - z 0 4- eo) + 2 fi(u - v’y - w’z) > (B7) 


a z = w + y0 


By substitution of the expressions (Al) for y and z, the accelerations 
of the particle are obtained as follows in terms of the coordinates of the 
cross-section £ and rj: 


= u - £2 2 (x + u) - 2 Q.v + tj [v’ cos p - w ' sin p + 

ft 2 (v' cos J 3 + w' sin p) + 2 J 20 sin p] + s[v f sin p - 
w’ cos p + Q 2 ( -v ' sin p + w’ cos p) + 2ft0 cos pj 

(v + e Q ) + 2ftu + rjj^-0 sin p - fl 2 (cos p - 0 sin 3) 
2n(-v’ cos p - w’ sin 3 0 + e[-$ COS P - 
i^C-sin p - 0 cos p) + 2Q(v’ sin 3 - w’ cos p)J 


a-^ = v - ft 


(B8) 


a z = w +■ -q0 cos p - £0 sin p 
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The inertia loadings on the heam may now he derived from these 
acceleration equations by appropriate integrations over the cross section. 
The total or resultant loadings desired in the analysis given in the body 
of this paper are the sum of the inertia loadings and the applied loadings, 
and are given by the following equations, when the assumption is made that 
the cross section is symmetrical about the major principal axis: 




where p is the density of the structural material and may be a function 
of £ and t). If the integrations are performed with the use of equa- 
tions (B8) , the following loadings are found; 
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p x = -m - & 2 (x + u) - 2ftv-j - me jj-v'cos 3 - w'sin 3 + 
n 2 (v'cos p 4- w T sin 3) + 2f20 sin p| 

Py = L y - mjy - ft 2 (v + e Q ) + - me |^0 sin 3 - 2 2 (cos 3 

0 sin 3) 4- 2ft(-v'cos 3 - v'sin p)j 

• • 

p z = L z - m(w) - me 0 cos 3 


q! x = M - |(v + e Q )sin 3 + e o 0 cos pj + me(v sin 3 - 

w cos 3 4- 2ftu sin 3 ) - j^k^ 2 - kmi 2 ) sin 3 cos 3 4- 

( k m 2 2 " W 2 )^ cos 2 ^| ~ 0 “^ m 2 " 2 to ( k m 2 2 “ ^ml 2 ) v ‘ cos £ sin 

2 fftn(k mP 2 sin 2 3 4- 2 cos 2 3 )w t 

q y = -Q^e [(x + u)sin 3 4- x0 cos pj 4- me(u sin 3 - 

20,-v sin p) + j^k^p 2 - k ml 2 )v , sin 3 cos 3 4- 4- 

l%l 2 cos 2 3)w|j - rn^kmp 2 - Wi 2 ) v'sin 3 cos 3 4- 


(BIO) 

(BH) 

(B12) 

3 - 
(B15) 


m(kmp 2 Bin 2 3 + l%ii 2 cos 2 3) (2n<z5 - w’) 


(B14) 
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q z = -ft 2 me [(x + u)cos 3 - x0 sin 3] + me(u cos 3 - 2ftv cos 3) + 

A "(k^ 2 - k^^w'sin 3 cos 3 + ^k m2 2 cos 2 3 + Tc ml 2 sin 2 3) v'J - 
m ( k m 2 2 " 3 cos 3 “ mfk TT ,p 2 cos 2 3 + k ml 2 sin 2 3^v' +- 

2to ( k m2 2 - W 2 )^ sin 3 cos 3 (B 15 ) 


Equations (BIO) to (BI5) give the general expressions for the 
loadings and contain many terms of second order which, for most engineering 
purposes can be neglected.. In equation (B 1 Q), for example, all other 
terms are small in comparison to mft 2 x, which is the conventional expres- 
sion for p x '(note that p x ± -H?'); for most - practical purposes these 
small terms may therefore be neglected. For any specific application, 
however, particularly in the case of unusual configurations, the relative 
magnitude and importance of the terms should be determined. In the case 
of conventional helicopter or propeller blades, it_is believed that the 
following first-order reductions of equations (BIO) to (BI5) are suitable 
for most engineering applications : - 


p x = -T' = m£l 2 x 

Py = Ly - mjlr - Q 2 (v + e Q )J + me Jp sin 3 + & 2 (cos 3-0 sin I^J 
p z = L z - m^w + e*0* cos 3 ^ 
q x = M - fi^e £(v +- e 0 )sin 3 +-e o 0 cos pj + 

me(v sin p - w cos 3 ) - f&n ^k mP 2 - kml 2 )(sin 3 cos 3 + 

0 cos 23) - m 

q = -Q 2 mex(sin 3 -K0""cos p) 
q = -fi 2 mex(cos P - 0 sin 3 ) 


(B 16 ) 


where 
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(b) Rotating beams. 

Figure 1.- Case’s treated in previous theory. 
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^ y-Centroid 

Elastic axis 
(a) Coordinates. 



(b) Displacement 



(c) Moments. 

Figure 3«- Nomenclature. 
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Figure 7-- Natural frequencies of example 2 showing the importance of 

centrifugal-force coupling. 


